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Eddy  Activity  in  the  California  Current 
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abstract 

Two  cABairiAlc  tourca  for  the  feoc.'Atioe  o'  maotc*}*  eddy  Mlinty  ic  '.Et  CaJiforXLi*  Current  «;«  iocaI 
IwxaeUlue  iaetobUity  Aod/or  (be  wLad  t’.ro*  Ad>iAeat  lo  tbe  co«*t.  TVe  Utu:  resele 

witb  eddic*  propafatlsi  watwoed  from  tb«  co«<t  isto  ibe  CedrforruA  Currest  n*  Kottoy  <r»rt  dye.Ar&.ct 
la  Ihii  liudy,  two  viad  dhvee  modeit  Arc  utr'ised  to  '.eat  sbt  rtiAti»c  •ifatScAOct  of  ihoM  two  •o-ortco,  ooc 
AB  eddy  raolvisf  quAti*  gcoetropbic  (Q-C)  model,  wbid:i  bo*  tbe  AbUity  to  icprewest  bArocLaic  iaeiibvhly 
but  aot  tbe  eoaitAl  rtepooM  to  wiadi,  Aad  tbe  otbe:  a  i-l/T  Uyer  prirtuttire  eijvAtioB  {PEj  mode;,  wb.cb 
ba<  tbe  Ability  to  npreatst  tbe  cOAJtAl  reapccic  to  w-.ade  but  aot  bArtxIiaie  iaetAbiUty.  Both  model*  Lav* 
tbe  AAtnt  epAtiAl  (rid  (i.c.,  AppreximAUty  20  e.m)  ead  Are  drivea  by  tbe  aasx  courM-fnd  wAd  etrcM  fore^sf 
field*  ovei  tbe  lAme  oae  ytAi  time  period  (;  e..  Kovembe;  1818  lo  October  HtT),  tbre  period  ciceec  oecAiAe 
of  tbe  AVAliAbility  of  CtoeAt  Altimetrie  ica  level  obeervAtioa*  with  wbtcb  to  verify  tbcee  model*  EArLer. 
Wbite  ef  al.  (1880)  AaAlyicd  iboe  lAmt  eltimetric  ma  level  obeervAtloa*.  fiadlai  domlaeat  m*»oACAle  eddy 
Activity  oecurrisi  oo  wAveleaftb  kaIc*  of  e0&-800  km  Add  period  Acalet  of  6-12  mosib*.  propA(Atiat  to  tbe 
wot  At  2-S  em/(«e,  faitet  At  lower  lAtltvide,  loveroed  by  Ro**by  wove  dy&A.*alt»  We  fisd  ia  ibi*  *tady 
tbAt  tbe  PE  model  ii  able  lo  limuUte  ^uiliutively  ibi*  dittributlos  of  tbe  eddy  vAriAact  a*  it  api^cat*  ta 
Altimetrie  tea  level,  yleldlo(  litaificASt  cobereace  And  pbeje  betwrea  model  Asd  obaerred  wa  level  rtsiduAi* 
Aloni  loa(itudc/time  mAtricei  At  30* S  aad  40*  A*.  Tbe  Q-G  model,  oa  the  otber  bend,  it  foaad  iacApobic 
of  limulAting  even  tbe  quAlitAtive  dUtributloa  of  eddy  vAriADce.  Tbe  reoAOD  for  tbe  Afrecmeat  betweea  tbe 
PE  model  end  tbe  tAteilite  Altimetrie  cca  level  obeexvAtioai  i*  Cbet  tbe  deminAat  wrurce  of  mcAOACAle  rody 
Activity  os  tbeae  time  Aod  epAcc  icAle*  i*  tbe  wiad  foreisg  edyAcent  to  tbe  cooet,  nsoeitfitd  by  both  Roeeby 
And  Kelvin  wive  dysAmic*. 


1.  Introduction 

Mcsoscalc  eddy  activity  in  the  California  Current 
system  has  been  examined  repcaiedly  over  ihc  years 
(e.g.,  Wyllie  1966:  Bernstein  et  al.  1977;  Simpson  ct 
ai.  1984;  Rcinecker  et  al.  1987)  using  both  in  siiu  and 
satellite  observations.  The  study  of  eddy  activity  in  ihc 
California  Current  began  with  the  analysis  of  in  situ 
hydrographic  data  collected  during  a  sequence  of 
CALCOn  cn  ises  beginning  in  1950.  with  eddies  and 
fronts  embedded  in  the  California  Current  repeatedly 
observed  over  a  period  of  many  years.  In  long-term 
monthly  mean  dynamic  height  maps  from  CaLCOFI, 
Wyllie  ( 1966)  demonstrated  the  existence  of  semiper¬ 
manent  eddies  in  the  California  Current  off  Point 
Concepiion  and  adjacent  to  Punta  Eugenia.  Hickey 
( 1979)  later  demonstrated  the  existence  of  a  semiper¬ 
manent  eddy  west  of  Monierey  in  these  same  data. 
Prior  10  satellite  observations,  iransicni  eddy  aciivity 
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was  difficult  to  studs  wiih  ihe  C4LCOFI  hvdrograph 
observations  because  of  lack  of  adequate  lirric  st 
quences  ihat  would  allow-  individual  eddies  lo  b 
tracked  over  space-time.  The  advent  of  satellite  ob 
scrvaiions  made  possible  a  S'noptic  \  ie*  of  the  surfaci 
eddy  field,  often  repeated  for  weeks  and  months  at  : 
lime. 

Transient  eddy  activity  in  the  California  Currcrii  was 
first  studied  using  AVHRR  observations  by  Bernstein 
el  ai.  ( 1977).  finding  transient  mesoscale  eddy  activiiy 
in  the  California  Current  to  have  onginalcd  in  ncar- 
coastal  waters.  This  observation  was  followed  bs  Owen 
( 1980)  who  examined  mcsoscalc  eddy  actisity  In  the 
California  Current  from  a  sanety  of  in  situ  and  saiellue 
observations,  finding  semipermanent  eddies  associated 
with  coastal  and  bathymetric  irregularities  (eg.  the 
semipermanent  eddy  off  Point  Conception ),  and  tran¬ 
sient  eddies  associated  with  wmd  C'cnis  adjacent  to 
the  coast.  A  thorough  hydrographic  study  of  one  spe¬ 
cific  transient  eddy  observed  in  saicllilc  data  was  con¬ 
ducted  by  Simpson  ct  al.  ( 1984).  in  which  its  origin 
was  examined  in  the  light  of  theories  of  topographic 
generation  and  baroclinic  instability.  This  eddy  was 
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found  to  have  coastal  undercurrent  water  at  its  center, 
suggesting  its  origin  at  the  coast.  These  satellite  studies 
began  to  suggest  that  transient  mesoscale  edd>  activav 
in  the  California  Current  had  us  origin  m  the  near 
coastal  region,  propagating  westward  from  there  into 
the  current. 

A  sequence  of  localized  intensive  hvdrographic  sur¬ 
veys  in  a  small  subregion  of  the  California  Current 
region  was  conducted  (called  OPTOM,<K  )  m  order  to 
study  the  behavior  of  mesoscale  eddv  activity  in  the 
current  (e.g.,  Robinson  et  al.  1986;  Rcmecker  ei  al. 
1987).  The  hi^-resoiuiion  data  of  the  OPTOM.^  ex¬ 
periment  was  inserted  by  Reinecker  et  al.  ( 198’)  into 
a  quasigeostrophic  numerical  model,  performing  dy¬ 
namical  interpolation  and  providing  forecasts  of  the 
eddy  evolution  in  the  survey  region.  They  found  these 
mesoscale  eddies  tracking  westward  through  the  subre¬ 
gion  at  Rossby  wave  speeds.  Based  on  an  analysts  of 
the  energy  budget,  they  found  that  forcing  by  the  curl 
of  the  wind  stress,  as  well  as  baroclinic  instability  pro¬ 
cesses.  were  essential  for  understanding  the  evolution 
of  the  observed  eddies  in  the  current.  This  work  had 
little  to  say  about  the  origin  of  these  eddies,  focusing 
rather  on  their  modification  by  the  local  wmd-sircss 
cur!  and  the  background  current. 

Mesoscale  eddy  activity  over  the  entire  California 
Current  region  was  examined  by  White  et  al.  ( 1990). 
based  upon  altimetric  sea-level  observations  from  the 
Geosat  Exact  Repeat  Mission  (ERM)  over  the  one- 
year  period  from  January  1987  to  December  1987. 
Analysis  of  sea-level  residuals  from  the  long-term  mean 
permitted  a  visualization  of  the  space-time  evolution 
of  the  mesoscale  eddy  activity  over  the  entire  region 
not  previously  possible.  The  rms  differences  of  the  al¬ 
timetric  sea  level  about  the  mean  showed  maximum 
eddy  activity  at  three  principal  locations  adjacent  to 
the  coast  of  California;  that  is.  off  Puma  Eugenia  at 
ZT’N.  southwest  of  Point  Conception  at  32'’N.  and 
adjacent  to  the  coast  between  Monterey  and  Cape 
Mendocino.  Time-longitude  plots  demonstrated 
coastal  mesoscale  eddy  activity  propagating  westward 
into  and  through  the  California  Current. 

Numerical  models  arc  powerful  tools  in  the  inves¬ 
tigation  of  mesoscale  eddy  generation  and  dynamics. 
A  good  e,xample  of  this  is  given  by  Holland  and  Lin 
(I975a.b)and  Holland  ( 1978),  who  used  both  quasi¬ 
geostrophic  (QG)  and  primitive  equation  ( PE)  models 
in  studying  the  role  of  mesoscale  eddies  in  the  general 
circulation  of  the  ocean.  Holland  established  that  me¬ 
soscale  eddies  spontaneously  arise  from  instabilities  in 
the  mean  ocean  currents  and  are  subsequently  con¬ 
trolled  by  nonlinear  Rossby  wave  dynamics.  Recently, 
numerical  models  have  developed  the  capability  for 
conducting  realistic  simulations  that  allow  this  and 
other  hypotheses  to  be  tested.  A  particularly  relevant 
modeling  study  of  the  variability  of  the  eastern  Pacific 
is  the  one  done  by  Cummins  et  al.  ( 1 986 )  using  a  sim¬ 
plified  one-layer  QG  model.  They  analyzed  the  simu¬ 
lated  Rossby  wave  field  in  the  eastern  Pacific  and  found 


It  to  be  dominaicd  b>  wa'.ei  emanating  hv'm  the  east¬ 
ern  boundary,  organized  into  iVihciciit  patterns  I  las 
showed  >cr>  eicafiv  ihal  variable  v^ir-.d  .vifi  ncai  ibc 
coast  can  be  a  wiurce  tor  Rossbv  lor  ibe  whoic 

eastern  Pacific  The  model  used  m  itos  'anr.f  studs  dij^rs 
not  haxe  the  capacity  otticcuming  bofs>ciiiin.  iv  unstaMe 
(being  a  one-laser  mode! }  and  hence,  touid  not  assess 
the  importance  of  this  contnbution  to  ihc  oscraii  eddy 
energy  budget  of  the  region  More  recent  siudies  teg, 
Auad  Cl  al.  199!  S  base  shown  that  for  the  QG  model 
in  genera!  the  contnbution  so  ihe  eddy  sanabihty  in 
the  California  Cuncni  region  from  the  wind-euri  forc¬ 
ing  IS  overshadowed  by  the  consnbuiion  from  baro- 
chnic  instability 

Therefore,  tsxo  sources  for  the  genfraiion  of  meso- 
Kalc  eddy  activity  m  the  California  Curreni  arc  hy¬ 
pothesized.  I  c..  ItKai  barcKlintc  insiabihty  of  the  main 
current  and  remote  generation  of  eddy  activity  by  ihe 
wind  stress  adjacent  to  the  coast,  the  latter  [vropagaung 
westward  into  land  through!  the  California  Current 
via  Rossby  wave  dynamics  To  test  these  hypotherers, 
two  numencal  mtvdcis  arc  utilized  One  nn.xic!  is  eight 
layers,  eddy  rcsol'  mg  and  quasigeostrophic  ( ()C> ).  w  uh 
the  ability  to  represent  harcKlimc  msiabiiiiy  hut  not 
the  coastal  response  to  winds  The  other  model  is  a 
1  i/;-laycr  primitive  equation  (PEL  with  the  ability  to 
represent  the  coastal  response  to  winds  but  not  baro¬ 
clinic  instability.  Both  numencal  mcxlels  have  a  similar 
spatial  gnd  ( i.e..  approximately  20  kmj  and  arc  drnen 
by  the  same  coarse-gnd  wind  Stress  forcing  fields  over 
the  same  one-year  time  penod  (i  c  .  January  1987  to 
December  1987).  Each  mode!  aiicmpis  to  provide  a 
realistic  simulation  of  the  mesoscale  eddy  activity  m 
the  California  Current.  These  model  results  arc  then 
compared  with  those  observed  during  the  last  year  oi 
the  Geosat  ER.M  (White  et  al.  1990).  tc.siing  the  eddy 
generation  mechanism  isolated  m  each  model.  Of 
course,  each  of  these  models  ( C?G  and  PE)  include 
"more  physics"  than  just  baroclinic  instability  and  local 
wind  forcing  at  the  coast.  However,  as  argued  below, 
for  the  California  Current  area  and  for  the  realistic 
wind  used,  those  two  mechanisms  dominate  the  re¬ 
sponse  of  the  QG  and  PE  model,  respectively. 

2.  Geosat  altimetric  sea-lc*el  residuals  in  the 

California  current 

Geosat  altimetric  sea-level  residuals  m  the  California 
Current  region  for  the  period  January  1 987  to  Decem¬ 
ber  1987  arc  used  m  this  study  in  a  comparison  with 
the  model-generated  mesoscale  eddy  activity,  .A  com¬ 
plete  description  of  the  Geosat  Exact  Repeat  Mission 
(ERM)  in  this  region  is  given  by  White  ct  al.  (  1990). 
in  this  study,  only  a  brief  account  of  the  ERM  and  the 
extensive  prcprcKCSsing  of  the  raw  data  is  given.  This 
is  followed  by  a  discussion  of  the  rms  of  the  inira- 
annual  variability  about  the  1987  mean,  taken  from 
White  et  a!.  ( 1990). 

The  Geosat  ERM  has  a  17-day  repeat  orbit  cycle 
and  a  longitudinal  track  separation  of  approximately 
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i.47  degrees.  Along  ihe  ground  irack.  one-second  av¬ 
erages  of  altimeiric  sea-level  measuremenis  arc  formed 
with  a  spatial  resolution  of  approximately  6  7  km  The 
grid  resulting  from  the  superposition  of  all  the  ascend¬ 
ing  and  descending  orbits  m  the  California  Current 
region  is  shown  in  Fig.  1.  A  one-year  penod  (from 
January  1987  to  Oecember  1987)  is  available  for  use 
in  this  study,  previously  analyzed  by  White  et  al. 
(1990).  Preprocessing  of  the  altimctnc  sea-level  ob¬ 
servations  begins  with  an  application  of  the  cnviron- 
mcnul  corrections  supplied  with  the  raw  observations. 
Next,  to  remove  the  gcoid  and  residual  orbit  error,  the 
track  mean  is  subtracted  from  each  individual  obser¬ 
vation,  producing  altimeiric  sca-leve!  residuals,  and  the 
resulting  residuals  are  detrended  along  each  track.  A 
Gaussian  filter  is  then  applied  to  the  alongirack  alti- 
mctric  sea-level  residuals,  suppressing  the  vanance  of 
wavelengths  less  than  50  km.  This  latter  procedure  is 
conducted  as  a  quality  control  c.xercise  to  suppress 
white  noise  assumed  to  be  of  instrumental  origin.  Fi¬ 
nally.  the  alongtrack  data  is  decimated  at  25-km  in¬ 
tervals.  Only  the  ftltered/decimatcd  altimctnc  sea-level 
residuals  along  ascending  tracks  are  used  in  this  study. 

These  ftltered/decimatcd  altimetric  sea-level  resid¬ 
uals  are  mapped  onto  a  regular  0.5®  latitude-longitude 
grid  (White  ct  al.  1990)  over  the  California  Current 
region  every  17  days  over  the  one-year  period  of  in¬ 
terest.  The  subsequent  spatial  distribution  of  rms  dif¬ 
ferences  about  the  one-year  mean  is  displayed  in  Fig, 
2.  Overall,  maximum  rms  variability  is  confined  to  the 
coastal  ocean  in  three  different  latitude  locations.  A 
latitude  maximum  occurs  adjacent  to  the  coast  of 
Northern  California  from  37°N  to  40“ N.  extending 
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Fig.  I.  Distribution  of  the  ascending  and  descending  tracks  from 
the  GEOSaT  Exact  Repeat  Mission.  Each  track  was  separated  Trom 
adjacent  tracks  by  approximately  140  km;  each  track  was  repeated 
every  1 7  days  for  the  one-year  period  (January  1987-December  IW) 
of  this  study. 


Fig  ’  Observed  spaua)  distisbuucn  ot  ihe  rmr  rcssduah  c! 
altimctrw  lea  k'Cl  aboui  ihe  ont.^car  mean  from 
December  i9S7 


from  .Monterey  lo  Cape  Mendocino  Another  Icxa! 
maximum  occurs  southwest  of  Point  Conception  near 
32“N.  120'W''.  where  frequent  mcsoscalc  eddy  activity 
(i.e.,  the  Southern  California  Eddy)  was  noted  rcpcat- 
cdly(c.g..  Simpson  ct  al.  1986  ).  A  third  local  maximum 
occurs  d’rectly  west  of  Punta  Eugenia  where  frequent 
mcsoscalc  eddy  activity  was  also  noted  earlier  (eg. 
Wyllie  1966).  Most  importantly  wuh  regards  to  this 
study  are  the  extensions  of  these  rms  variability  max¬ 
ima  toward  the  west,  suggesting  westward  propagation 
of  mcsoscalc  eddy  energy  from  these  apparent  source 
regions,  each  extending  from  ihc  coast  out  into  the 
California  Current  itself 


3.  Primitive  equation  model  of  the  California  current 

The  U/;  layer  PE  mcxici  used  in  this  study  is  de¬ 
scribed  in  detail  by  Pares-Sierm  and  O  Brien  ( 1989) 
It  consists  of  one  dynamically  active  layer  of  constant 
density  p  and  variable  depth  //  over  an  infinitely  deep 
lower  layer  of  higher  density  p  Ap.  The  model  do¬ 
main  extends  from  20"-50“N  and  from  the  coast  of 
North  America  to  150“W,  The  eastern  boundary  of 
the  model  domain  realistically  follows  the  geometry  of 
the  eastern  North  Pacific  coastline.  The  nonhern. 
southern,  and  western  boundaries  arc  open,  vvnh  a 
Sommerfeld  radiation  condition  ; c.g..  Camcrlcngo  and 
O'Brien  1980)  imoosed  there.  The  spatial  resolution 
over  this  model  domain  is  approximately  20  km  in 
both  the  zonal  and  meridional  directions.  The  equa¬ 
tions  of  the  model  arc 
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where  8  and  4>  are  ihc  latitude  and  longitude,  respec¬ 
tively;  l  and  Fare  the  transports  in  the  east-west  and 
nonh-south  directions,  respectively;  H  is  the  depth  of 
the  upper  layer;  g'  is  reduced  gravity;  r*  and  t*  arc  the 
wind-stress  components;  ,4  is  the  eddy  viscosity  coef¬ 
ficient;  and  a  is  the  radius  of  the  eanh.  The  values  of 
the  parameters  are  given  in  appendix  A. 

Earlier,  Pares-Sierra  and  O'Brien  ( 1989 )  had  forced 
this  nniodel  with  realistic  monthly  mean  wind  stress  (i  c,. 
computed  from  the  COaDS  surface  wind  observation 
set),  to  simulate  the  interannuai  variability  of  the  Cal¬ 
ifornia  Current,  The  pattern  of  the  average  circulation 
in  the  Noah  Pacific  derived  from  the  1 1/:  layer  PE 
model  was  demonstrated  to  be  similar  to  that  observed. 
Moreover,  upper-layer  seasonal  and  interannuai  vari¬ 
ability  in  the  model  was  demonstrated  to  be  similar  to 
that  observed  at  a  number  of  coastal  sea  level  stations 
{ Pares-Sierra  and  O’Brien  1989),  At  that  time,  Pares- 
Sierra  and  O’Brien  ( 1989)  found  most  of  the  seasonal 
variability  in  the  model  originating  at  or  near  the  coast. 
Subsequent  off-shore  propagation  of  this  coastal  vari¬ 
ability  appeared  to  provide  for  significant  mesoscaie 
eddy  aaivity  in  the  California  Current  { Pares-Sierra 
1991 ).  Thissupponed  the  earlier  suggestion  by  Mysak 
( 1983)  that  annual  fluctuations  at  the  eastern  boundary 
were  responsible  for  most  of  the  mesoscaie  eddy  activity 
in  the  interior  ocean. 

In  the  present  study,  the  Fleet  .Numerical  Oceano¬ 
graphic  Center  (FNOC)  Surface  Wind  Analyses  arc 
used  to  compute  synoptic  wind-stress  estimates  for 
driving  the  PE  model  over  the  period  coincident  with 
the  Geosat  ERM.  The  synoptic  FNOC  wind-stress  es¬ 
timates  have  a  resolution  of  two  degrees  in  the  merid¬ 
ional  and  zonal  directions,  and  of  one  day  in  lime; 
these  estimates  arc  intenxjiated  down  onto  the  10  km- 
20  minute  space-time  grid,  upon  which  the  model  is 
integrated.  The  model  integration  is  conducted  in  three 
stages:  the  model  is  driven  from  rest  for  lO  years  with 
the  annua!  mean  wind-stress  estimates;  then,  the  model 
is  driven  for  another  ten  years  with  the  long-term  mean 
annual  cycle  of  the  wind  stress;  finally,  the  model  is 
driven  for  ten  years  from  1978  to  1987  with  the  syn¬ 
optic  wind-stress  estimates.  .Model  results  from  the  final 
year  (i.e..  January  l987~December  1987)  arc  com¬ 
pared  with  the  Geosat  altimeiric  sea-level  observations. 

Earlier,  Pares-Sierra  and  O'Brien  (1989)  demon¬ 
strated  that  the  large-scale  features  of  the  California 
Current  at  the  sea  surface  can  be  simulated  in  quah- 
tative  fashkin  with  a  wind-driven  reduced  gravity  1  i/i 
layer  PE  model.  This  is  demonstrated  in  the  upper 
panel  of  Fig.  3,  where  the  mean  upper-layer  thickness 
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presence  of  a  coastal  countercurreni.  sornetimes  called 
the  Davidson  Current.  The  Davidson  Current  is  stron¬ 
gest  in  winter  along  the  northern  portion  of  the  Cali¬ 
fornia  coast.  The  California  Current  is  seen  to  be  well 
separated  from  the  coast  of  California;  hence,  (he  two 
sources  of  mesoscale  eddy  activity  to  be  e.xamined  in 
this  study  (i.c..  baroclinic  insiabiliiy  and  coastal  wind 
generation )  arc  also  seen  to  be  separated  geographically. 

In  the  PE  model,  the  location  of  maximum  vari¬ 
ability  in  response  to  synoptic  FNOC  wind-stress  forc¬ 
ing  is  not  in  the  California  Current  itself  but  m  the 
waveguide  adjacent  to  the  coast  of  California.  This  can 
be  seen  in  the  lower  panel  of  Fig.  3.  w  here  the  spatial 
distribution  of  the  rms  of  sea-level  residuals  is  dis¬ 
played.  Overall,  the  maximum  rms  vanability  is  con¬ 
fined  to  the  coastal  ocean  in  two  different  latitude  lo¬ 
cations.  A  latitude  maximum  occurs  adjacent  to  the 
coast  of  Northern  California  from  37*  to  40*.S'.  ex¬ 
tending  from  .Vlontcrcy  to  Cape  Mendocino,  very  sim¬ 
ilar  to  that  observed  in  Fig.  2.  .Another  local  maximum 
occurs  southwest  of  Point  Conception,  where  intense 
variability  is  observed  in  the  local  wind-stress  field  east 
of  the  Channel  Islands.  These  areas  of  enhanced  me¬ 
soscale  variability  coincide  with  those  found  by  Cum¬ 
mins  ct  al.  { 1986).  Most  important  arc  the  extensions 
of  the  two  rms  variability  maxima  in  Fig.  3  toward  the 
southwest,  consistent  with  that  observed  in  Fig.  2  and. 
again,  suggestive  of  westward  propagation  from  the 
coast  into  the  California  Current  located  in  the  offshore 
region. 

This  rms  distribution  is  the  result  of  an  annual  pat¬ 
tern  of  upwelling  and  downwelling  associated  with  the 
annual  north-south  reversal  of  the  wind  stress  along 
the  coast.  Superimposed  on  this  arc  the  upwelling/ 
downwelling  events  associated  with  synoptic  changes 
in  w’ind-stress  forcing;  these  synoptic  events,  however, 
are  generally  much  smaller  than  those  associated  wiih 
the  annual  cycle  of  wind-stress  variability.  The  latter 
results  from  a  deepening  (weakening)  of  the  semiper¬ 
manent  continental  thermal  low  over  California  dunng 
spring-summer  (autumn-winter),  which  generates 
strong  northwesterly  (weak  northwesterly )  wind  stress 
parallel  to  the  coast.  Strong  northwesterly  wind  stress 
is  upwelling  favorable,  represented  in  the  PE  model  by 
a  reduction  in  the  depth  of  the  upper  layer.  Weak 
northwesterly  wind  stress  along  the  southern  pan  of 
the  domain,  and  a  reversal  of  the  direction  of  the  wind 
stress  in  the  vicinity  of  Cape  Mendocino,  during  au¬ 
tumn-winter  weakens  the  upwelling  pattern.  This  is 
interspersed  with  regions  of  coastal  downwelling.  Sea¬ 
sonal  coastal  variability  is  strongest  in  the  vicinity  of 
Cape  Mendocino,  w-here  an  actual  seasonal  wind-stress 
reversal  is  present,  with  a  secondary  maximum  occur¬ 
ring  southwest  of  Point  Conception. 

4.  Eddy-resolving  quasigeostrophtc  model  of  the 

California  current 

A  quasigeostrophic  (QG)  eddy-resolving  model 
(Holland  and  Vallis  1990)  is  used  to  investigate  the 


role  of  baroclinic  instability  as  a  ntechanssm  for  gen¬ 
erating  mesoscale  eddy  acto-ay  m  the  Cahlorrua  Cur¬ 
rent.  Details  ol  the  phySics  ofihis  can  lx;  found 

in  Holland  (  197g).  and  discussion  of  ssmilamies  and 
differences  between  molnlaycr  QG  and  PE  models  is 
conducted  in  Holland  (  1985.  1986)  The  sa.mc  mixJcl 
is  used  by  Auad  et  al  {  1991 )  to  study  the  circulation 
and  energetics  of  the  eastern  Pacific  from  the  point  of 
view  of  QO  dynamics.  The  equations  of  the  model  arc 
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is  the  (5G  poicni’sl 'omciiy  for  the  5  !h  layer  ( i.  »  l. 
8 ).  Here  V  x  t  is  the  w md-stress  curl,  p,,  is  the  reference 
density  of  the  fluid.  //;  is  the  i  ih  layer  thickness.  6, 
is  a  Kronccker  delta.  An  is  the  lateral  fnction  cocfficieni. 

is  the  total  sireamfunclion  for  ihc  A  ih  layer,  t  ts  the 
bottom  fnction  coefficient.  C  and  arc  the  Conolis 
parameter  and  its  meridional  gradient,  respect ociy. 
taken  at  the  model  central  latitude  r,  =  33'N.  h{x.  r) 
is  the  bottom  topography.  The  total  dcnvatisc  operator 
is  DlDt  =  d/di  +  {  )).  where  the  last  term  is 

the  Jacobian  operator.  The  fu.,  :  symbol  represents 
the  venicai  displacement  of  the  (A  f';!ih  interface, 
which  can  be  regarded  as  an  isopycnal  surface; 
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where  gl.,  ;  is  the  reduced  gravity  corresponding  to 
the  (A  -t-  i/;)th  interface.  An  impoaant  characteristic 
of  this  model,  in  regard  to  the  present  study,  is  the 
form  that  the  boundary  condition  takes  m  terms  of  the 
primary  field  The  boundary  condition  imposed  at 
the  lateral  boundaries  is  that  of  no  normal  How.  that 
is. 

/j .  f  =  0  1^  =  const  al  the  boundary; 

in  particular,  the  depth  of  the  upper  layer  (h  x  4^, 
-  4'-})  is  constant  along  the  boundary  (hut  a  function 
oftimc)(  Holland  1978).  The  method  for  solving  these 
equations  has  been  given  in  many  papers  (e  g..  Holland 
1978;  Cummins  and  Mysak  1988;  Holland  and  Vallis 
1990). 

This  QG  model  extends  over  the  same  domain  as 
the  PE  model  (i.e..  20®-50*N.  coast  ofN’orth  .Amcnca- 
I50®W).  with  approximately  the  same  spacc-timc  grid 
resolution  and  with  eight  layers  in  the  vertical. 

Many  earlier  studies  using  this  QG  model  concen¬ 
trate  upon  simulating  mesoscale  eddy  activity  m  the 
very  energetic  western  boundary  currents  (e  g.,  &hmitz 
and  Holland  1986).  In  the  study  of  the  California  Cur¬ 
rent.  a  different  set  of  model  parameters  is  chosen,  re¬ 
quiring  the  model  to  simulate  caslcrn  boundary  current 
variability  with  much  greater  vcnical  resolution  in  the 
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upper  500  m  of  ocean  and  wiih  less  density  ct  >  ;rasi 
between  layers.  The  mode!  parameters  used  m  this 
study  arc  given  in  appendix  B.  The  reduced  gravity 
estimates  are  based  on  historical  CalCOFI  hydrographic 
observations  taken  near  the  California  coast. 

Integration  of  this  eddy-resolving  QG  model  occurs 
in  step-wise  fashion  similar  to  that  used  in  integrating 
the  PE  model.  The  mean  wind-stress  curl  used  to  force 
the  model  is  computed  from  the  FNCXT  Surface  Wind 
Analysis  extending  for  nearly  10  years  from  1977  to 
1987.  The  QG  mode!  is  integrated  from  rest  using  the 
annual  cycle  of  monthly  mean  wind-stress  curl  until 
statistical  equilibrium  is  reached;  then  the  model  is 
driven  through  three  cycles  of  the  7-y  car  synoptic  wind- 
stress  curl.  During  the  last  year  of  this  integration  ( Jan- 
uary-Dccember  1987),  QG  model  upper-layer  thick¬ 
ness  estimates 


i.e.,  h  =  -r-  -  4':) 

are  saved  for  a  comparison  between  those  observed 
and  those  from  the  PE  model. 

The  spatial  distribution  of  mean  upper-layer  depth, 
for  the  January- December  period  from  the  eddy-re¬ 
solving  QG  model  is  given  in  the  upper  panel  of  Fig. 
4.  TTiis  distribution  is  qualitatively  similar  to  that  given 
in  the  upper  panel  of  Fig.  3  obtained  from  the  1  i/y 
layer  PE  model,  since  both  maps  are  determined  by 
the  distribution  of  the  mean  wind-stress  curl.  In  each 
case,  the  California  Current  is  found  to  occur  250-750 
km  off  the  coast  of  California. 

The  spatial  distribution  of  the  rms  variability  in  sea 
level  from  the  eddy- resolving  QG  model  is  shown  in 
the  lower  panel  of  Fig.  4.  The  QG  model  presents  a 
very  different  spatial  pattern  of  variability  than  that  of 
the  PE  model  ( and  of  the  ohserved  data )  shown  in  the 
lower  panel  of  Fig.  3.  Most  conspicuous  is  the  absence 
of  strong  variability  along  the  coast,  particularly  be¬ 
tween  Monterey  and  Cape  Mendocino.  Rather,  the 
dominant  characteristic  of  the  QG  model  reflected  in 
the  spatial  distribution  of  rms  is  its  ability  to  sustain 
baroclinic  instability  in  the  mean  currents.  This  is  most 
strongly  evidenced  by  the  presence  of  a  zonal  band  of 
high  rms  variability  in  the  southern  portion  of  the 
model  domain  (20®-25‘'N).  This  occurs  in  the  region 
of  westward  return  flow  of  the  subtropical  gyre,  where 
conditions  for  baroclinic  instability  arc  more  favorable 
than,  say.  for  the  California  Current  (e.g.,  Lee  1988). 
Yet,  to  a  lesser  degree,  baroclinic  instability  occurs  in 
the  California  Current  as  well,  evidenced  by  an  along¬ 
shore  band  of  maximum  values  of  the  rms  differences 
superimposed  upon  the  mean  position  of  the  California 
Current.  Auad  et  al.  ( I99l ),  using  this  same  model. 
showed  the  highly  energetic  eddy  fields  of  the  coastal 
region  (their  area  I)  to  be  fed  by  energy  transmitted 
fi^m  the  wind  into  available  potential  energy  of  the 
mean  flow;  subsequently,  baroclinic  instability  pro¬ 
cesses  transform  available  potential  energy  into  eddy 
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kinetic  energy.  Their  analysis  also  showed  hjriHiinic 
instability  dominating  barolropic  insiahiiiucs  m  !his 
cddy-gcncraiing  process.  As  pointed  out  by  \ujd  ei  a! 
(1991),  these  energy  paths  represent  only  .oeraecd 
processes.  The  contention  that  baroclinic  msiability 
dominates  barotropic  instability  or  direct  wind  foremi; 
does  not  mean  that  these  other  processes  are  absent 
Other  sources  of  mesoscalc  variability  arc  oi  course 
possible  in  the  QG  formulation  (and  the  PH  nuHlel ) 
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besides  baroclinic  instability.  Two  candidates  are 
barotropic  instability  and  direct  wind-stress  curl  vari¬ 
ability  at  the  eastern  boundary.  In  the  study  of  Cum¬ 
mins  et  al.  ( 1986)  the  latter  mechanism  provides  ail 
the  variability  in  the  model.  In  their  work  waves  can 
clearly  be  seen  emanating  from  the  eastern  boundary 
al  two  principal  locations  near  the  coast  and  propa¬ 
gating  southwestwardly  toward  the  far  field.  Although 
this  mechanism  of  generation  is  certainly  present  in 
the  QG  model  used  here,  its  response  is  overshadowed 
by  the  effect  of  baroclinic  instability,  as  demonstrated 
by  Auad  ct  al.  ( 1 99 1 ) . 

Another  difference  between  the  PE  and  QG  model 
is  the  freedom  of  the  former  to  respond  to  the  irrota- 
tional  part  of  the  wina  stress  that  is  precluded  a  priori 
in  the  QG  model  By  decomposing  the  wind-stress  field 
into  an  irrotational  part  and  a  nondivergent  part  [i.e.. 
r  =  +  r' ,  where  V  x  *  0  and  V  •  r''  =  0,  sec 

Morse  and  Feshbach  (1953)],  we  can  quantify  the 
contribution  ( not  shown )  of  the  irrotational  part  of  the 
wind  in  the  PE  model.  By  forcing  the  PE  model  with 
only  t'  we  find  that  the  response  in  the  open  ocean  is 
minimal.  At  the  coast,  although  an  evident  seasonal 
cycle  exists,  its  magnitude  is  negligible  compared  to 
the  magnitude  of  the  seasonal  cycle  produced  by  the 
nondivergent  pan  of  the  wind. 

The  main  dynamical  difference  between  the  two 
models  is,  however,  the  filtering  of  gravity  waves  in  the 
QG  model.  As  discussed  below,  this  is  the  decisive  fac¬ 
tor  explaining  the  lack  of  observed  mesoscale  variability 
in  the  coastal  QG  model.  Moreover,  it  explains  the 
inability  of  the  QG  model  to  simulate  the  observed 
radiation  of  coastal  eddy  energy  into  the  ocean  interior. 
3  process  that  wc  argue  is  determinant  in  the  eastern 
Pacific. 

S.  Intercomparison  between  observed  and  model 

RMS  variability 

The  charaaeristic  spatial  distribution  of  the  rms  re¬ 
siduals  in  altimciric  sea  level  displayed  in  Fig.  2  are 
compared  with  those  of  the  PE  and  QG  models  dis¬ 
played  in  Figs.  3b  and  4b.  As  discussed  earlier,  maxi¬ 
mum  rms  variability  in  observed  altimetric  sea  level 
occurs  adjacent  to  the  coast  of  Nonhem  California 
from  37'’N  to  40®N.  coincident  to  the  region  of  max¬ 
imum  annual  variability  of  the  wind  stress  (c.g..  Nelson 
1977).  Southwest  of  Point  Conception,  a  local  maxi¬ 
mum  in  rms  variability  exists;  a  local  minimum  in  the 
rms  occurs  in  the  California  Bight;  and  another  local 
maximum  occurs  adjacent  to  Punta  Eugenia. 

To  affect  the  comparison  of  the  spatial  distribution 
of  the  two  sets,  the  sea-level  residuals  from  the  PE 
model  are  treated  to  exactly  the  same  preprocessing 
filters  as  applied  to  the  altimetric  sea-level  observations 
described  in  section  2.  Also,  the  two  datasets  are  taken 
over  exactly  the  same  time  period.  Upon  inspection  of 
these  two  spatial  distributions,  that  from  the  PE  "^odel 
is  seen  to  be  qualitatively  similar  to  that  observed,  t  he 


rms  spatial  distnbution  of  the  PE  model  d!vpia>s  a 
maximum  between  Montcres  and  Cape  Mendoesno. 
similar  to  that  observed  m  both  configuration  and 
magnitude.  Farther  to  the  south,  the  PE  model  repro¬ 
duces  the  local  rms  ma.ximum  observeo  southwest  of 
Point  Conception.  Again,  a  local  rms  minimum  in  the 
California  Bight  separates  the  former  maximum  from 
the  coast. 

The  contours  of  rms  variability  in  the  PE  model  are 
concentrated  somewhat  closer  to  the  coast  than  wuh 
the  observed  altimetric  sea  Level,  particularly  north  of 
35‘’N.  This  is  probably  due  to  the  suppression  of  west¬ 
ward  propagation  of  coastal  disturbances  nonh  of  ap¬ 
proximately  35°N  on  period  scales  less  than  one  year, 
that  is.  to  the  existence  of  a  critical  latitude  ( McCreary 
et  al.  1987).  in  the  particular  PE  model  wc  are  using 
( 1 1/2  layers),  coastal  variability  is  not  allowed  to  prop¬ 
agate  into  the  California  Current  north  of  35®N.  except 
on  time  scales  exceeding  one  y  ear.  Mesoscale  eddy  ac¬ 
tivity  in  the  model  California  Current  north  of  35'N 
is  associated  with  periods  that  exceed  one  year  in  du¬ 
ration.  This  is  only  partially  true  in  the  obserx aliens, 
where  both  longer  and  shorter  period  activity  is  found 
in  the  California  Current  nonh  of  35°N.  A  finer  modal 
structure  in  the  model  is  probably  needed  to  simulate 
this  broadened  spectral  band. 

The  spatial  distnbution  of  the  rms  sca-lcvei  residuals 
from  the  cddy-resoKing  QG  model  is  shown  in  Fig. 
4b.  Upon  inspection  of  observed  and  QG-model  spatial 
distributions  of  rms  (Fig.  2  and  4b).  the  QG  model 
can  be  seen  to  present  a  very  different  spatial  pattern 
of  variability  than  is  observed.  Most  conspicuously  the 
QG  mode!  is  unable  to  simulate  the  intense  variability 
along  the  coast,  particularly  between  Monterey  and 
Cape  Mendocino.  Rather,  the  spatial  distribution  of 
rms  variability  in  the  QG  model  favors  the  generation 
of  variability  in  the  major  currents  that  constitute  the 
subtropical  gyre  (i.e..  the  California  Current  and  the 
.North  Equatorial  Current).  This  is  induced  by  baro¬ 
clinic  instability  processes.  This  pattern  of  variability 
in  the  observed  rms  vanabiliiy  is  overshadowed  by 
coastal  seasonal  variability.  Moreover,  the  magnitude 
of  these  rms  sea-level  residuals  in  the  California  Current 
are  smaller  than  those  observed,  suggesting  that  local 
baroclinic  instability  processes  in  the  California  Current 
arc  dominated  by  mesoscale  eddy  activity  generated  at 
the  coast. 

6.  Coherence  and  phase  between  the  observed  and 

model  sea  level 

Having  shown  that  the  PE  model  i.s  able  to  .stmulaic 
the  spatial  distribution  of  the  rms  sca-lcvcl  residuals 
about  the  mean,  wc  now  examine  the  model's  ability 
to  simulate  the  phase  of  the  dominant  observed  sca- 
icvcl  variability.  Longitude-lime  matrices  for  observed 
and  model  sea-level  residuals  are  prc.scnicd  at  40° N  in 
Fig.  5  and  at  30°N  in  Fig.  6.  Each  longitude-time  ma¬ 
trix  extends  from  ihc  coast  of  California  westward  to 
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140®W  for  the  one-year  period  January  1987-Deceni- 
bcr  1987.  The  spectral  coherence  and  phase  between 
the  two  matrices  at  each  latitude  arc  also  given. 

At  40°N,  both  the  observed  and  model  sea-level  re¬ 
siduals  are  dominated  by  a  strong  annual  cycle  (upper 
panels  Fig.  5).  From  November  to  May.  values  arc 
positive  along  the  coast,  indicating  a  downwelling  re¬ 
gime  brought  on  by  strong  poleward  wind  stress.  After 
May,  positive  values  propagate  westward  vrith  negative 
values  developing  at  the  coast.  In  both  matrices,  the 
region  of  positive  and  negative  residuals  propagates 
westward  from  the  coast  as  the  residuals  of  opposite 
sign  develop  at  the  coast.  This  process  is  associated 
with  Rossby  wave  propagation  ( Pares-Sierra  1991). 
Note,  however,  in  both  matrices  this  westward  prop¬ 
agation  is  arrested  at  250-500  km  from  the  coast,  as¬ 
sociated  with  a  rapid  reduction  in  the  magnitude  in 
the  offshore  direction.  In  the  PE  model,  this  is  a  man¬ 
ifestation  of  the  Rossby  wave  critical  latitude,  where 
waves  of  annual  period  generated  at  the  coast  are 
trapped  near  the  coast  at  40®N.  this  due  to  their  strong 
offshore  decaying  factor  at  these  periods  ( McCreary  et 
al.  1987).  From  the  good  qualitative  comparison  be¬ 
tween  these  two  longitude-time  matrices  in  the  coastal 
waveguide,  we  can  infer  that  similar  dynamical  pro¬ 
cesses  are  operating  in  the  observed  situation  as  well 
to  trap  eddy  energy  in  the  coastal  waveguide.  West  of 
the  coastal  waveguide  ( i.e,.  in  the  California  Current ). 
both  the  observed  and  model  longitude-time  matrices 
display  longer  period-scale  variability  required  to  es¬ 
cape  the  coastal  trapping  mechanism.  The  observed 
matrix  does  contain  higher-period  mesoscale  eddy  ac¬ 
tivity  west  of  the  coastal  waveguide,  but  its  energy  is 
small  compared  to  the  longer  period  activity  that  has 
its  origins  at  the  coast;  its  source  could  be  baroclinic 
instability  or  local  forcing  by  the  time-varying  wind 
curl. 

Within  the  longitude-time  matrices  at  40°N,  ob¬ 
served  mesoscale  eddy  activity  displays  a  w-ider  range 
of  spatial  and  temporal  scales  than  in  the  PE  model. 
Model  sea-level  residuals  are  heavily  dominated  by  the 
annual  frequency,  with  little  influence  of  the  higher- 
frequency  variability  that  is  observed.  We  believe  this 
to  be  due  to  the  lack  of  baroclinic  instability  in  the 
model,  possibly  exacerbated  by  the  lack  of  mesoscale 
synoptic  wind-stress  estimates  with  which  to  force  the 
PE  model.  The  differences  in  the  patterns  of  spectral 
energy  density  in  the  zonal  wavenumber-frequency 
domain  ( not  shown )  confirm  this  inability  of  the  model 
to  simulate  the  wider  bandwidth  of  spectral  response 
to  the  wind-stress  forcing.  For  the  model  sea-level  re¬ 
siduals.  the  spectral  energy  density  is  concentrated  at 
frequencies  lower  than  about  0.3  cycles/month  (i.e.. 
corresponding  to  periods  greater  than  approximately 
4  months).  For  the  observed  sea-level  residuals,  most 
of  the  spectral  energy  is  also  below  0.3  cycles/ month, 
but  above  this  frequency  the  spectrum  is  relatively  flat, 
with  an  order  of  magnitude  greater  spectral  energy 
density  than  in  the  PE  model. 


Both  zonal  wavenumber-frequency  spectra  at  40“N 
display  peaks  in  their  spectral  energy  density  in  the 
negative  wavenumber  quadrant  near  the  theoretical 
Rossby  wave  dispersion  curve,  the  latter  calculated  us¬ 
ing  the  I  i/2-taycr  PE  model  characteristics.  This  indi¬ 
cates  that  most  of  the  energy  in  both  the  real  and  model 
oceans  propagates  westward  in  the  form  of  qua  si-nnear 
Rossby  waves,  with  shorter  periods  adjacent  to  the  coast 
and  longer  penods  farther  offshore. 

Squared  spectral  coherence  between  the  observed 
and  model  sea-level  residual  matrices  at  40°N  (bottom 
panel  of  Fig.  5)  is  significant,  with  squared  coherence 
larger  than  0.5  for  the  range  0.0  to  0.03  cycies/day  and 
-0.5  to  -0.17  cycies/degree  of  longitude.  Therefore, 
at  40°N  the  PE  model  can  explain  more  than  70%  of 
the  variance  in  that  zonal  wavenumber-frequency 
range.  These  scales  of  vanabiliiy  also  dominate  the 
zonal  wavenumber-frequency  spectra  and  provide  a 
definition  in  this  study  for  mesoscale  eddy  activity. 
Spectral  phase  indicates  that  for  this  range  of  zonal 
wavenumbers  and  periods,  where  high  coherence  exists, 
the  model  and  observed  signals  are  nearly  in  phase. 

.At  30°N.  the  ability  of  the  wind-driven  PE  model 
to  simulate  the  observed  allimctric  sca-le'cl  residuals 
(top  panel  in  Fig.  6)  is  basically  the  same  as  at  40' N 
( Fig.  5 ).  Significant  spectral  coherence  occurs  between 
the  observed  and  model  sea-level  residuals  over  a  range 
of  frequencies  that  include  annual  and  semiannual  pe¬ 
riods.  with  a  phase  difference  of  less  than  45®  (i.e..  less 
than  1.5  months).  One  difference  is  that  at  these  period 
scales,  westward  propagation  of  signals  from  the  coast 
is  not  arrested  by  the  influence  of  the  cntical  latitude. 
Moreover,  in  agreement  with  linear  Rossby  wave  the¬ 
ory,  westward  propagation  of  these  signals  from  the 
coast  is  faster  at  30® N  than  at  40° N.  as  characterized 
by  a  larger  slope  of  the  positive  and  negative  contours 
in  both  longitude-lime  matrices  in  Fig.  5  compared  to 
those  in  Fig.  6. 

8.  Discussion  and  conclusions 

In  the  Introduction,  we  hypothesized  two  sources 
for  the  generation  of  mesoscale  eddy  activity  in  the 
Caliiomia  Current.  They  arc  local  baroclinic  instability 
of  the  current  and/or  remote  generation  of  eddy  ac¬ 
tivity  by  the  wind  stress  adjacent  to  the  coast,  the  latter 
propagating  westward  into  (and  through)  the  California 
Current  via  Rossby  wave  dynamics.  In  this  study,  two 
models  are  utilized  to  test  the  relative  significance  of 
these  two  sources.  One  model  is  eddy  resolving  and 
quasigeostrophic  (QG).  with  the  ability  to  represent 
baroclinic  Instability  but  not  the  coastal  response  to 
winds,  "^hc  other  model  is  a  noneddy  resolving  1 1/;- 
laycr  prii...tive  equation  ( PE),  with  ihc  ability  to  rep¬ 
resent  the  coastal  response  to  winds  but  not  baroclinic 
instability.  Both  models  have  similar  spatial  grids  (i.e.. 
approximately  20  km)  and  are  driven  by  the  same 
coarse-grid  wind-stress  forcing  fields  over  the  same  one- 
year  Time  period  (i.e.,  January  1987  to  December 
1987).  This  period  is  chosen  because  ofthe  avaifabilily 
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of  Geosat  altimeiric  sea-level  obsenaiions  with  which 
to  verify  the  models. 

We  find  in  this  study  that  the  PE  model  is  able  to 
simulate  qualitatively  the  observ'cd  distribution  of  the 
eddy  variance  as  it  appears  in  altimetric  sea  level.  It 
also  yields  significant  coherence  and  phase  between 
modeled  and  observed  sea-level  residuals  along  longi¬ 
tude-time  matrices  at  30°N  and  40°N.  The  QG  model, 
on  the  other  hand,  is  found  incapable  of  simulating 
the  dominant  qualitative  features  of  the  eddy  variance. 
The  reason  for  the  agreement  between  the  PE  model 
and  the  Geosat  altimetric  set-level  observations  is  that 
the  dominant  source  of  mesoscale  eddy  activity  on 
these  time  and  space  scales  is  the  wind  forcing  adjacent 
to  the  coast.  Therefore,  we  conclude  that  the  wind- 
driven  coastal  generation  of  mesoscale  eddy  activity 
dominates  that  of  baroclinic  instability  in  the  produc¬ 
tion  of  mesoscale  eddy  activity  in  the  California  Cur¬ 
rent. 

The  more  significant  difference  between  the  QG 
model  and  the  PE  model  in  relation  to  this  study  is 
the  lack  of  the  necessary  degrees  of  freedom  of  the 
former  to  sustain  gravity  waves  [i.e..  Poinare  and  Kel¬ 
vin  waves.  Ripa  ( 1991 )].  Although  Poincare  waves  are 
probably  not  relevant  for  the  scales  being  addressed 
here,  the  filtering  of  Kelvin  waves  [by  both  the  as¬ 
sumption  of  subineniality  (i.e.,  1  /cr/ 1 ),  as  well  as 
by  the  requirement  of  constancy  of  (and  of  /i)  along 
the  coast]  probably  suppresses  a  large  percentage  of 
coastally  generated  variability  in  the  model.  At 
least,  the  contribution  to  the  mesoscale  variability  gen¬ 
erated  equatorward  of  some  area  of  interest  would  be 
absent.  This  latter  source  has  been  shown  to  be  very 
important  at  certain  scales  (e.g.,  Pares-Sierra  and 
O’Brien  1989).  Even  more  importantly,  the  QG  re¬ 
quirement  of  constant  layer  depth  along  the  coast  very 
effectively  buffers  the  response  of  the  coastal  ocean  to 
localized  wind  events.  Wind  forcing  would  have  to  be 
very  strong  and  the  events  would  have  to  last  long 
enough  to  lift  the  whole  boundary  at  the  same  time. 

We  recognize  that  QG  models  are  adequate  for 
studying  some  processes  in  the  California  Current,  but 
one  of  these  processes  is  not  the  generation  of  the  total 
mesoscale  variability.  One  of  its  key  ingredients, 
namely  direct  wind  forcing  of  coastal  variability  is 
missing  from  the  physics  available  in  this  model.  It  has 
been  shown  in  this  study  that  the  eastern  Pacific  Ocean 
is  not  a  suitable  location  for  utilization  of  the  QG  model 
physics.  On  the  other  hand,  if  we  were  trying  to  study 
the  dynamics  of  baroclinic  instability  in  the  California 
Current,  for  example,  the  opposite  would  be  true;  that 
is,  the  1 1/2-layer  PE  model  would  be  useless  and  the 
QG  model  would  be  very  appropriate. 
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Summary  of  Rc<fucod-Gra>ity  Vtodcl  Parami-icrs 


Coetficicnt  of  eddy  mscoshv 
Drag  coefficient 

Reduced  gravity 

Initial  upper-layer  depth 

350  m 

1.5  V  itr^ 

0,01  m  s'- 
200  m 

APPENDIX  B 

Summary  of  Quasigeost. . 

vhic  .Viodel  Parameters 

Layer  thicknesses  (m) 

100.  150.  250.  450.  650.  900. 

1 150.  1350 

Reduced  gravities  (m  s'') 

0.007061.  0008196,  0.008389. 
0.005521. 0.002  PO. 
0.0005094,  0.00006908 

Lateral  friction  coefficient 
(m-  s") 

200 

.fo(s-‘) 

7.943  X  10'* 

#f(m*'  s'") 

1.92  X  lO'” 

Earth's  radius  (m) 

6.37  X  10* 
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